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Abstract 



We report a search for the flavor-changing neutral current decay B — ► 
using a 29.1 fb 1 data sample accumulated at the T(45) resonance 
with the Belle detector at the KEKB e + e~ storage ring. We observe the decay 
process B — > K£ + £~(£ = e,fi), for the first time, with a branching fraction of 
B(B -► K£+£~) = (0.75l°;lf ± 0.09) x 10" 6 . 

PACS numbers: 11.30.Hv, 13.20.He, 13.25.Hw 
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Flavor-changing neutral current (FCNC) processes are forbidden at the tree level in the 
Standard Model (SM), but are induced by loop or box diagrams. If non-SM particles partic- 
ipate in the loop or box diagrams, their amplitudes may interfere with the SM amplitudes. 
This makes FCNC processes an ideal place to search for new physics. 

The b — > s transition is a penguin-diagram mediated FCNC process. The CLEO group 
reported the first observation of the B — > X s 7 radiative penguin decay jl). The measured 
branching fraction for this process has been used to set the most stringent indirect limit on 
the charged Higgs mass and to constrain the magnitude of the effective Wilson coefficient 
of the electromagnetic penguin operator \Cj S \ 0. However, it cannot constrain the sign 
of Cj S , which is essential to obtain definitive evidence of new physics since C| ff is negative 
in the SM while it can be positive in some non-SM physics models The electroweak 
penguin decays B — ► X s £ + £~~ are promising from this point of view since the coefficients 
Cg and Cxo can be determined by measuring the dilepton invariant mass distributions 
and forward-backward charge asymmetry of the dilepton and the B — > X s ^y decay rate 

Standard Model branching fraction predictions for B — > K^*'£ + £~ decays are listed in 
Table | Although several groups || have searched for B — > K^*H + £~ decays, no 

evidence has been observed. 

In this Letter, we present the results of a search for B decays to K^*H + £~ using data 
collected with the Belle detector at the KEKB storage ring || . The data sample corresponds 
to 29.1 fb~ x taken at the T(45) resonance and contains 31.3 million BB pairs. 

Belle is a general-purpose detector based on a 1.5 T superconducting solenoid magnet 
that surrounds the KEKB beam crossing point. Charged particle tracking is provided by a 
Silicon Vertex Detector (SVD) and a Central Drift Chamber (CDC). Particle identification 
is accomplished by a combination of silica Aerogel Cerenkov Counters (ACC), a Time of 
Flight counter system (TOF) and specific ionization measurements (dE/dx) in the CDC. A 
CsI(Tl) Electromagnetic Calorimeter (ECL) is located inside the solenoid coil. The h/Kl 
detector (KLM) is located outside of the coil. A detailed description of the Belle detector 



can be found elsewhere [ 10] . 



In this analysis, charged tracks, except for the K s — > 7t + tt~ decay daughters, are required 
to have a point of closest approach to the interaction point within 0.5 cm in the r(p plane 
and 5.0 cm in the z direction, where the r<p plane is the plane perpendicular to the electron- 
beam (z) direction. Electrons are identified from the ratio of shower energy in the ECL 
to the momentum measured by the CDC, the shower shape of the cluster in the ECL, 
dE/dx in the CDC and the light yield in the ACC. Tracks are identified as muons based 
on the matching quality and penetration depth of associated hits in the KLM. To reduce 
the misidentification of hadrons as leptons, we require that the momentum be greater than 
0.5 GeV/c and 1.0 GeV/c for electron and muon candidates, respectively. Charged kaons 
and pions are identified by a likelihood ratio based on dE/dx in the CDC, time-of-flight 
information and the ACC response. 

Photons are selected from isolated showers in the ECL with energy greater than 50 MeV 
and a shape that is consistent with an electromagnetic shower. Neutral pion candidates 
are reconstructed from pairs of photons, and are required to have an invariant mass within 
10 MeV/c 2 of the nominal 7r° mass and a laboratory momentum greater than 0.1 GeV/c. K s 
candidates are reconstructed from oppositely charged tracks with a vertex displaced from 
the interaction point. We require the invariant mass to lie within 15 MeV/c 2 of the nominal 
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Kg mass. 

K* candidates are formed by combining a kaon and a pion: K + tt~, K^n , Kg7i + or K + tt 
TL\ . The K* invariant mass is required to lie within 75 MeV/c 2 of the nominal K* mass. 
For modes involving 7r°'s, combinatorial backgrounds are reduced by the further requirement 
cos^hd < 0.8, where 6*h e i is defined as the angle between the K* momentum direction and 
the kaon momentum direction in the K* rest frame. 

B candidates are reconstructed from a K^*> candidate and an oppositely charged lepton 
pair. Backgrounds from the B — > J/ip(ip are rejected using the dilepton invariant mass 
veto windows; -0.25 < M ee - M J/4 , < 0.07 GeV/c 2 for J/ip K*, -0.20 < M ee - M J/ ^ ) < 
0.07 GeV/c 2 for J/ij, K$ K&), -0.15 < M w - Mj^ < 0.08 GeV/c 2 for J/ip K* and 
-0.10 < M m - Mj/^i) < 0.08 GeV/c 2 for J/ip K{$ To suppress the background 

from photon conversions and tt° Dalitz decays, we require the dielectron mass to satisfy 
M ee > 0.14 GeV/c 2 . 

Backgrounds from continuum qq events are suppressed using event shape variables. A 



Fisher discriminant T [0 is calculated from the energy flow in 9 cones along the B candidate 



sphericity axis and the normalized second Fox- Wolfram moment R2 |13| . Furthermore, we 
use the B meson flight direction cos 8b and the angle between the B meson sphericity axis 
and the z axis, cos# sp h. For the muon mode, cos# sp h is not used since its distribution is 
nearly the same for signal and continuum due to detector acceptance. We combine J 7 , cos 6b 
and cos^ph into one likelihood ratio CR cont defined as £TZ cont = £ sig /(£ si g + C coa t), where 
£ sig and £ con t ar e the products of the probability density functions for signal and continuum 
background, respectively. 

The major background from BB events is due to semileptonic B decays. The missing 
energy of the event, E miss , is used to suppress this background since we expect a large amount 
of missing energy due to the undetected neutrino. The B meson flight angle cos 9b is also 
used to suppress combinatorial background in BB events. We combine E miss and cos 9b into 
the likelihood ratio £TZ BB -, defined similarly to ClZ cont . 



Finally, we calculate the beam-energy constrained mass M bc = J -E 2 eam — p% and the 
energy difference AE = E B — E hcam to select B candidates, where E heam = yfs/2 is the 
beam energy in the center of mass (cm) frame and ps and Eb are the measured momentum 
and energy of the B candidate in the cm frame, respectively. The selection criteria are tuned 
to maximize the expected significance S/ V S + B where S is the signal yield and B is the 
expected background in the signal box. S and B are determined from GEANT based Monte 
Carlo (MC) samples. The B -> K^*H + t decays are generated according to the Greub, 
Ioannissian and Wyler model with the branching fractions predicted by Ali et al. 
The interference between K^*H + £~ and J/ip(ip is not considered. The signal box is 

defined as |M bc — M B \ < 7 MeV/c 2 (2.7c) for both the electron mode and the muon mode, 
where Mb is the nominal B meson mass, and —0.06 < AE < 0.04 GeV for the electron 
mode and |Ai?| < 0.04 GeV for the muon mode. We make selections on ClZ corit and CR-bb 
that reject 85% of the continuum background and 45% of the BB background and retain 
75% of the signal for all modes except for those with Kgn + and K + n° final states, where 
the selection on CR, b ~b ^ s tightened to reject 55% of the BB background and retain 70% of 
the signal. The overall detection efficiencies, estimated by the MC simulation, are listed in 
Table |TJ. 

To determine the signal yield, we perform a binned maximum-likelihood fit to each 
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Mb c distribution. The expected number of events is calculated as a function of Mb c , from 
a Gaussian signal distribution plus background functions. The mean and the width of 
the signal Gaussian are determined using observed J/ipK^ events. A MC study shows 
that the width has no dependence on the dilepton invariant mass. The background from 
real leptons is parameterized by the ARGUS function flT4" |. The shape is determined from 



400 fb _1 MC samples, each containing at least one oppositely charged lepton pair. As 
shown in Figure [l| (right column), the ARGUS function is a good representation of the 
background distributions. The MC shape is consistent with the shapes derived from the 
AE sideband and the K^e ± fx ¥ samples in the data. The background contribution due 
to misidentification of hadrons as muons is parameterized by another ARGUS function 
and a Gaussian. The ARGUS function represents the combinatorial background while the 
Gaussian represents the background that makes a peak in the signal box. The shape and 
normalization of this background are fixed using the B — > Kh + h~ data sample (h refers to 
hadrons). All Kh + h~ combinations are weighted by the momentum dependent probability 
of misidentifying Kh + h~ as K{x + fi~ . This study yields 0.27 ± 0.03 Kh + h~ events in the 
peak region. For electron mode, the misidentified Kh + h~ background in the peak region 
is less than 0.007 events. Other backgrounds with misidentified leptons are negligible. The 
normalizations of the signal and the background from real leptons are floated in the fit. 

The fit results are shown in Figure [I] (left column) and summarized in Table The 
statistical significance is defined as ^—2 ln(£ /£max), where £ max is the maximum likelihood 
in the M^ c fit and Co is the likelihood when the signal yield is constrained to be zero. We 
observe 11 Kfx + /j~ events. The fit to the Mb c distribution yields 9.5i3^ signal and 1.6 ±0.4 
background events. The statistical significance of this excess is 4.7. The probability of an 
upward fluctuation of the background to 11 or more events is 5.5 x 10~ 6 , which corresponds 
to 4.4 standard deviations for a Gaussian probability distribution. As a test we also perform 
a fit to the AE distribution and find a signal yield of 8.5i2!L which is consistent with the 
Mb c fit results. 

The kinematical properties of the Kfi + fi~ events are further examined to check for po- 
tential backgrounds that might peak in the signal area. The B + — > D°tt + , D° — > K + tt~ 
decay chain is the largest expected source of Kh + h~ background. We expect 0.20 ± 0.12 
events based on a MC simulation study. Another possible background source is double- 
misidentification of the B — > J/tpK, J ftp —>■ fi + fi~ decay chain where the kaon and a 
muon are misidentified as a muon and a kaon, respectively. The K + fi~ combinations 
with K + tt~ and yU + /i~ hypotheses are examined for the candidate events, and show no 
cluster in the D° mass or J/ip mass region, which confirms the MC expectation. The 
B — > J/ipX, J/ip — > decay chain can be another background source when muon pairs 

from J/i[) — > fi + fi~ decays evade the J/ip veto. We expect 0.08 events using a MC sample. 
The fi pair effective mass distribution is consistent with the MC expectation (Fig. |^), and 
we observe no events close to the J/ip or if>' veto region. To summarize, we observe no 
indication of a background producing a peak in the M^ c distribution in the Kfi + fi~ sample. 

We consider systematic effects from the fit and the efficiency determination. Uncertainty 
in the background function is the dominant source of the systematic error. To evaluate 
the effect of the signal function parameters, the mean and the width of the Gaussian are 
changed by ±lcr from the values determined from J/ipK^*' events. The uncertainty in the 
background shape is obtained by varying the ARGUS shape parameter by ±lcr from the 
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value determined with a large MC sample. The magnitude of the variation is rescaled to an 
equivalent luminosity of 29.1 fb _1 . Even if the background shape is modified to maximize 
the background contribution in the signal region, the statistical significance of the K[i + 
signal remains above 4.0. The systematic errors associated with the fit function are shown 
in the third column of Table [TT]. Systematic uncertainties on the tracking, charged kaon 
ID, charged pion ID, electron ID, muon ID, K s detection and tt° detection efficiencies are 
estimated to be 2.3 to 2.5%, 2.1 to 2.5%, 0.8%, 1.8%, 2.2%, 8.7% and 6.8% per particle, 
respectively. 

In calculating the branching fraction, we assume equal fractions of charged and neutral 
B meson pair production at the T(AS). We combine neutral and charged I?-meson results 
for B — > Kfi + [i~ modes and obtain the branching fraction 

B(B -> Kfvr) = (0.99lgiHii) x 10- 6 , 

where the first and second errors are statistical and systematic, respectively. If we combine 
the B — > Ke + e~ and B — > Kji + \i~ decay modes, we observe 13.6l 3 J signal events with 
a statistical significance of 5.3. Assuming lepton universality, the branching fraction is 
determined to be 

B(B -> K£+r) = (0.75l°f ± 0.09) x lO" 6 . 

These values are consistent with the SM predictions For the modes with significance 

of less than 3.0, we also set upper limits for the branching fractions, employing the approach 
of Feldman and Cousins |15[], as listed in Table [TI|. These limits are consistent with SM 



predictions. 

In summary, we have observed the electroweak penguin decay B — ► K£ + £~. The branch- 
ing fractions obtained can be used to constrain contributions of new physics in the Wilson 
coefficients Cg S and C w . 
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ence, and Technology of Japan and the Japan Society for the Promotion of Science; the 
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sources; the Department of Science and Technology of India; the BK21 program of the 
Ministry of Education of Korea, the Basic Science program of the Korea Research Founda- 
tion, and the Center for High Energy Physics sponsored by the KOSEF; the Polish State 
Committee for Scientific Research under contract No.2P03B 17017; the Ministry of Science 
and Technology of Russian Federation; the National Science Council and the Ministry of 
Education of Taiwan; and the U.S. Department of Energy. 
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TABLE I. Branching fractions for B — > K^£ + £ predicted in the framework of the Standard 
Model. 



Mode 



Predicted branching fraction [ x 10 



Ali et al. [§ Greub et al. g| Melikhov et al. 

Kl+£- 0.57l£|{j 0.33 ± 0.07 0.42 ± 0.09 

K*e+e~ 2.3lg;| 1.4 ±0.3 1.4 ±0.5 

K*/j, + h~ l-9±o:i 1.0 ±0.2 1.0 ±0.4 



TABLE II. Summary of the fit results and branching fractions. Number of events observed 
in the signal box, number of signal and background events estimated from the M^ c fit, detec- 
tion efficiency of each mode, branching fraction obtained, 90% confidence level upper limit of the 
branching fraction and the statistical significance of the signal. The first error in the signal yield 
and branching fraction is statistical and the second one is systematic. The error in the efficiency 
includes MC statistics and systematic error. The error in the background is statistical only. 



Mode 


Observed 


Signal 


Back- 


Efficiency 


B 


U.L. 


Stat. 


events 


yield 


ground 


[%] 


[xl0~ 6 ] 


[xl0~ 6 ] 


signif. 


K°e + e- 


1 


n rr + 1. 4+0.4 
u -°-0.5-0.5 


0.3 ±0.2 


5.5 ±0.6 




2.7 




K + e + e- 
Ke + e~ 


5 
6 


o c+2.5+0.5 
8-0.7 
A -,+2.7+0.6 
^• 1 -2.1-0.8 


1.5 ±0.4 
1.7 ±0.4 


21.6 ±2.0 
13.6 ±1.3 


n r -,+0.37+0.09 
u - ol -0.27-0.11 
n ,,0+0.32+0.09 
u -^°-0.24-0.11 


1.4 
1.3 


2.4 
2.5 


K* e + e~ 


9 


A n+2.9 + 1.0 
^• u -2. 2-1.1 


3.9 ±0.7 


6.6 ±0.7 




6.4 




K* + e + e~ 
K*e+e~ 


4 
13 


9 C-+2.3+0.3 
^• d -l. 6-0.4 

C O+3.7 + 1.0 

°-°-3.0-l.l 


1.5 ±0.5 
5.7 ±0.9 


3.1 ±0.4 
4.8 ±0.5 


9 ns +1.23+0.35 
z - uo -l. 00-0.37 


8.9 
5.6 


2.5 


K+fi+ijr 


2 
9 
11 


i q+1. 8+0.0 
i - y -l.l-0.1 
7 O+3.4+0.9 
'■ -2. 7-1.0 
n cr+3. 8+0.8 
y -°-3. 1-1.0 


0.1 ±0.1 
1.8 ±0.4 
1.6 ±0.4 


6.5 ±0.7 
23.9 ±2.2 
15.2 ± 1.4 


n QA +0.88 +0.11 
u - y4 -0. 54-0.12 
fl qq+0.46+0.15 
u - yo -0. 36-0.16 
n QQ+0.40+0.13 
u - yy -0.32-0.14 


3.3 


2.8 
3.9 
4.7 




6 


o O+2.6+0.6 
°- z -l. 9-0.7 


2.2 ±0.5 


8.3 ±0.9 




4.2 




K*+fi+fi~ 
K*fj,+(j,~ 


2 

8 


n n +0.7+0.0 

u - u -o.o-o.o 

9 -1+2.9+0.9 
z - 1 -2.1-1.0 


2.7 ±0.6 
4.9 ±0.8 


3.5 ±0.4 
5.9 ±0.7 




3.9 
3.1 




K£ + £~ 


17 


io r+45+0.9 
i0 - D -3.8-l.l 


3.3 ±0.5 


14.4 ± 1.4 


°- 75 -a2i ± °- 0y 




5.3 
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FIG. 1. M bc distributions with fits for (a,b) B -> Ke + e~, (c,d) 5 -> K*e + e~, 
(e,f) 1? — »■ Kfi + fi~, (g,h) i? — »■ K*n + n~ and (i,j) Z? — »■ K£ + £~. Left column is for data and 
right column is for MC background. The solid curve in the right column shows the fit results with 
the ARGUS function, which is used in the fit to the data in the left column. 
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FIG. 2. Dimuon mass distribution of B — > K/j, + fj, candidates. The hatched histogram shows 
the data distribution while the open histogram shows the MC signal distribution. 
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